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1. Introduction 

The problem of Pomeron interaction in zero transverse dimensions have been discussed in the framework of the 
Reggeon calculus M about three decades ago. However, recently, we have seen a revival of the interest to this 
problem (see j^, ||, [| fj], ||, ^| and references therein). The very reason for this in our opinion is related to 
a hope to solve the old problem of finding the high energy asymptotic behaviour of the scattering amplitude in 
QCD . We hope for a solution not in the mean field approximation, where the solution has been discussed and well 
understood both analytically and numerically (see jnj [ll], |l^, [l3[ |lj, [ll| [l6| 0), but in the approach where the 



so called Pomeron loops should be taken into account |18, [L9L EG, Ell The problem of taking into account the 



Pomeron loops can be reduced to the BFKL Pomeron calculus |23j, E4l, t25| or/and to the solution of statistical physics 
problem: Langevin equation and directed percolation p6| , [27], |28| ] . The last approach is based on the probabilistic 
interpretation of the Pomeron calculus which also has roots in the past |2{| ^] . 

The BFKL Pomeron calculus in zero transverse dimensions being an oversimplified model has the same description 
in terms of the directed percolation as the general approach. Thus solving this model we can gain an experience that 
will be useful for the solution to a general problem of interaction of the BFKL Pomerons in QCD. 

It is well known that the BFKL calculus in zero transverse dimensions can be treated as a system that evolves 
in imaginary time it — Y with the Hamiltonian: 

H = —A(f><f) + + A (</>0 +2 — 4> 2 4> + ) and evolution equation for wave function — — = —H^ (1-1) 
v ' a Y 
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where the Pomeron intercept A oc as and the triple Pomeron vertex A oc a s . 

In the next section we will discuss the evolution equation for the generating functional that describes the system 
of Pomerons in terms of probabilities to find 'wee' partons (color dipoles [pl|). We introduce L(l — > 2) = A and 
r(2 — > 1) = A 7 with 7 is the amplitude for low energy interaction of the colour dipole with the target (target- 
Pomeron vertex) and 7 oc a s . The estimates for A and A are given in the leading order of perturbative QCD. We 



can trust the approach with the Hamiltonian of Eq. (1.1) only in the kinematic region of Y given by the following 
equation 

1 < ln(l/a|) < a s Y < — (1.2) 

a s 

Indeed, the n Pomeron exchanges give contribution which is proportional to (7 e Ay ) 2 and, therefore we need to sum 
them in the kinematic region where r ye AY > 1 or AY > ln(l/7) sa ln(l/a;|). However, we cannot go to ultra high 
energies since we do not know the higher order corrections to the BFKL kernel and to triple Pomeron vertex. The 
contribution of the BFKL Pomeron exchange can be written as ^y e ( A + Consta s) Y and the high order correction term 
is essential for as Y > l/as- 



In this paper for the kinematic range given by Eq. (1.2) we obtain two results. We introduce a new generating 



functional which allows us to calculate processes of multiparticle generation since it gives us the probability to have 
a given number cut and uncut Pomerons (see p2|). We derive the evolution equation for this generating functional 
both in the mean field approach and in the approach that takes into account the Pomeron loops. 

The second result is related to the method of summation of the Pomeron loops. We claim that in kinematic 



region of Eq. (1.2) the Pomeron interaction given by Hamiltonian of Eq. (1.1) can be reduce to the system of free 
Pomerons with the renormalized amplitude of dipolc-targct interaction at low energies. In other words, we can view 
the evolution of our system of 'wee' partons ( colour dipoles) as a system of not interacting partons only with emission 
absorbed in the evolution of the BFKL Pomerons, and all specific features of this system being determined by the 
low energy amplitude of 'wee' parton interaction with the target. Having this in mind we state that Mueller, Patel, 
Salam and Iancu |3^] approach gives the solution to the problem. 

The paper is organized as follows. In the next section we introduce the mean field approach and discuss it in 
the framework of the generating functional. We introduce a new generating functional which gives us a possibility 
to calculate the probability to find a given number of cut and uncut Pomerons. Therefore, knowing this generating 
functional we can calculate the cross section with given multiplicities. We derive the evolution equation for this 
functional. In section three we take into account the Pomeron loops and generalize the evolution equation. In this 



section we reconsider the problem of summation of all Pomeron loops in the kinematic region of Eq. ( L2 ) and argue 
that we can reduce this problem to consideration of a system of non-interacting Pomerons with renormalized vertex 
of Pomeron-target interaction. Based on this idea we use the Mueller, Patel, Salam and Iancu approach to calculate 
scattering amplitude at high energies both for elastic and inelastic interactions with different multiplicities of particles 
in the final state. 

In conclusions we summarize the results and discuss open problems. 

2. The mean field approximation 
2.1 General approach 

Our approach to multiparticle production is based on the AGK cutting rules pS^ . These rules stem from the 
unitarity constraint in s-channel, namely, 

2ImA BFKL (s,b) = 2N BFKL ( Sl b) = Gf n FKL {s,b) (2.1) 
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where ImA BFKL (s,b) = N BFKL (s,b) denotes the imaginary part of the elastic scattering amplitude for dipole- 
dipole interaction at energy W = yfs and at the impact parameter b. It is normalized in the way that the total cross 
section is equal to a tot — 2 J cPb N BFKL (s, b). Gf r FKL is the contribution of all inelastic processes for the BFKL 
Pomeron and <Ji n — J d b Gi 7l (s,b). Therefore, Eq. (2.1) gives us the structure of the BFKL Pomeron exchange 
through the inelastic processes and it can be formulated as the statement that the exchange of the BFKL Pomeron 
is related to the processes of multi-gluon production in a certain kinematics (see Fig. |l]). Eq. (2.1) is proven in p3[ 
in QCD. Using it, we can express all processes of multiparticle production in terms of exchange and interactions of 
the BFKL Pomerons and/or the cut BFKL Pomerons (see Fig. ||). 




a) b) c) 

Figure 2: Several examples of the Pomeron diagrams that contribute to the multiparticle production: the diffraction produc- 
tion of the bunch of particles in the region of rapidity ln(M 2 /m 2 ) = Y — Yo (Fig. []-a); the process of multiparticle production 
with the average multiplicity due to exchange of four Pomerons (Fig. ^-b); and the process of multiparticle production with 
the multiplicity in four times larger than the average multiplicity. 



The AGK cutting rules establish the relation between different 
processes that stems from BFKL Pomeron diagrams. For example, 
the simple triple Pomeron diagrams in Fig. |^ leads to three inelas- 
tic processes: the diffractive production of the system with mass 
ln(M 2 /m 2 ) = Y — Yo (Fig. || -A); the multi-gluon production in 
the entire kinematic region of rapidity Y — with the same multi- 
plicity of gluons as in one Pomeron ( Fig. ||-B); and the multi-gluon 
production in the region Y — but with the same multiplicity of 
gluons as in one Pomeron only in the rapidity window Y — Yo while 
for the rapidity Yq — the gluon multiplicity in two times larger 
than for one Pomeron( Fig. ||-C). The AGK cutting rules (32j say 
that the cross sections of these three processes are related as 




<ja ■ <?B 



1 



(2.2) 



At first sight the cross section of the process B is negative but it 
should be stressed that one Pomeron also contributes to the same 
process and the resulting contribution is always positive. 



Figure 1: The diagram for the cut BFKL 



Pomeron that illustrates Eq. (2.1) 
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Figure 3: The AGK cutting rules for the triple Pomeron diagram: the diffractive production of the system with mass 
ln(M 2 /m 2 ) = Y — Yo (Fig. ^ -A); the multi-gluon production in the entire kinematic region of rapidity Y — with the same 
multiplicity of gluons as in one Pomeron ( Fig. |^ -B); and the multi-gluon production in the region Y — but with the same 
multiplicity of gluons as in one Pomeron only in the rapidity window Y — Yo while for the rapidity Yo — the gluon multiplicity 
in two times larger than for one Pomeron( Fig. | -C) 



Fig. [j] and Fig. [| as well as Eq. (2.2) allow us to understand the equation for the single diffractive production in 
the mean field approximation that has been written by Kovchegov and Levin |34j and that has the following form 
(see Fig. ^ for graphical representation of this equation) 



B 




(^ D (x 02 ,b+ix 12 ,y,y ) + N D (x 12 ,b+^ Q 2,Y,Y ) - N D (x QU b,Y,Y ) 

+ iV c (x 02 ,b+ix 12 ,y,yo)^ I? (xi2,ixo2,y,lo) Eq. dU)(A) 

-4TV I3 (x 02 ,b+ix 12 ,y,y )iV (x 12 ,b+ix 02 ,y) Eq. Q(B) 



1 1 N 

+2jV (x 02 ,b + - Xl2 ,r)iVo(x 12 ,b+ -x 02 ,F) 



Eq. rt2j)(C) 



with the initial condition given by 



N D (x ± , b, Y = Y , Y ) = N$ (x x , b, Y ) 



(2.4) 



where Nq is the solution of the Balitsky - Kovchegov equation [|13| and Nd(xl, b; Y; Yo) is the diffraction dissociation 
of the colourless dipole with size x at impact parameter b into the system of gluon with the rapidity gap larger than 
Y at energy Y. At first sight, Eq. (2.3) contradicts the AGK relations given by Eq. fl2.2p (see Fig. pi), but this 



contradiction can be easily resolved if we take into account coefficient 2 in Eq. (2.1) (see |3jf for more details as well 
as for a proof based directly on the dipole approach to high energy scattering). 



Despite a simple structure of Eq. (2.3), which is only a little bit more complicated than the Balitsky - Kovchegov 
equation, as far as we know there exists no analytical solution to this equation and there is the only attempt to solve 
it numerically |$7|. However, this equation has a simple solution in the toy model (see |54| ^3, || ) which we are 
going to discuss. 



2.2 The BFKL Pomeron calculus in zero transverse dimensions: general approach 

The mean field approximation looks extremely simple in the BFKL Pomeron calculus in zero transverse dimensions 
(the toy model |3^, ^6|). Indeed, in the toy model, in which there is no dependence on the sizes of interacting 
dipoles, the generating functional degenerates to the generating function that has the form 



Z (u\Y) = Pn(Y)u n 



(2.5) 



n=0 



where P n {Y) is the probability to find n-dipoles (or/and n-Pomerons) at rapidity Y. For the probabilities P n (Y) we 
can write Markov chain, namely |3^, |3(J 



dP n (Y) 
dY 



-r(l -» 2) n P n (Y) + r(l -> 2) (n - 1) P„-i(Y) 



(2.6) 



Eq. (|j) has a simple structure: for the process of dipole splitting we see two terms. The first one with the negative 
sign describes a decrease of probability P n due to the process of splitting of dipoles. The second term with positive 
sign is responsible for the increase of the probability due to the same processes of dipole interactions. 



Eq. (2.6) can be re-written in the elegant form of the master equation for Zq, namely, 



dZ (u\Y) 



= -r(i ->2)u(i-«) 



dZ Q (u\Y) 



dY v ' v ' du 

where T(l — > 2) oc as in QCD. The initial and boundary conditions look as follows 



initial condition: 
boundary condition: 



Z (u\Y = 0) = u: 
Z (u = l\Y) = 1. 



(2.7) 



(2.8) 



In Eq. (2.8) the initial condition means that we are studying the evolution of one dipole, while the boundary condition 
follows from the normalization of the sum of probabilities. With this initial condition the linear differential equation 
Eq. (|2.7p can be written as non-linear one 



dZ (u\Y) 
dY 



= -r(l -> 2) Z (u\Y) + r(l - 2) Zl{u\Y) 



(2.9) 
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Introducing scattering amplitude 



N {~f\Y) = ImA el = -J2 



(-1)" d n Z (u\Y) j 
nl du ri 



I«=i7 



(2.10) 



with 7 being a scattering amplitude of interaction of single dipole with the target, we can find Eq. (2.7) for the 
amplitude |39| 

(2.11) 



^o(l\Y) =T{1 ^ 2){7 _^8N (j\Y) 



dY 



Using initial condition Nq{^\Y = 0) — 7 one can rewrite Eq. (2.11) as nondinear equation for the amplitude 

dN ( 7 \Y) 



dY 



r(l -> 2)7V ( 7 |y) _ r(l - 2)JV '( 7 |K) 



(2.12) 



Eq. (2.12) is easy to solve in this model and the solution has been found in |34], |38| |3|. It was noticed in |L3[ that 
one can get Eq. ( 2.12j ) directly from Eq. (2S) by substitution N (j\Y) = 1 — Zo(l — j\Y). Using this fact we can go 
back to Eq. (2/7) and write it as 



dY 



r(i-2) 7 



dZ(l—r\Y) 2 dZ(l- 7 \Y) 



(h. 



(2.13) 



2.3 The generating functional for the multiparticle production: definition and linear evolution equa- 
tion 

Here we would like to develop a method that will allow us not only to find the cross sections of the diffractive 
production but also to consider all processes of multi-particle production at high energy. Having this goal in mind, 
we propose a generalization of the generating functional given by Eq. (2.5), namely, we introduce a new generating 
functional Z(u,v\Y) as follows 

00 

Z{u,v\Y)= P{n,m\Y)u n v m (2.14) 

n— 0,m— 

where P (n, m\Y) is a probability to find n uncut Pomerons and m cut Pomerons. Directly from Eq. ( 2.14| ) and from 
the fact that P (n, to|Y) is a probability we find the first boundary condition 



Z(u= l,v= 1\Y) = 1 



(2.15) 



To find the second boundary condition we can use the full form of the s-channel unitarity constraint. Assuming 
that the scattering amplitude is pure imaginary at high energy, this constraint looks as follows 



2N(s,b) = \N(s,b)\ 2 + G v 



(2.16) 



where the first term in the l.h.s. is the elastic term with no cut Pomerons, while the second is the total contribution 
of the inelastic processes (in other words, sum over all cut Pomerons). Using functionals Zq(u\Y) and Z (u,v\Y) we 
can calculate the left and right hand sides of Eq. (2.16), namely, 



Ms, b)f 



N = 1 - Z a (l- 7 \Y) (see JT|) ; 
G m = 2 (1 - Z{1 - 7,1 - lm \Y)) 



(2.17) 
(2.18) 



where 7 = N(s = 0, b) is the imaginary part of the scattering amplitude of a dipole with the target at low 
energies, while ji n is the inelastic contribution (\N(s,b)\ 2 + Gt n ) for interaction of a dipole with the target at low 
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energy. Generally speaking both these amplitudes are arbitrary and have to be calculated from non-perturbative 
QCD, however, assuming the low energies are not very low and we can apply the relation of Eq. ( |2.l|) we see that 
2 7 = Jin- Using this relation we can re- write Eq. (2.17) and Eq. (2.18) in the form 



Z {u\Y) = Z(u,u\Y) 



(2.19) 



The initial condition depends on what we want to calculate. This is the main advantage of the generating 
functional that it allows us to calculate everything. For example, the cross section of single diffraction integrated 
over all produced masses (a s d) we can find just calculating Z (u,v = 0\Y — Yq) for the initial condition in the form 



Z(u,v\Y = Y ) 



The cross section is equal to 



o sd (j(Y )\Y-Y )) 



z(u = i-7,u = o|y-y ) 



(2.20) 



(2.21) 



The main idea of this paper is to introduce cut Pomeron being split into three different states. By analogy with 
Eq. ( [2.13 ) we can relate to each process correspondent term of differential equation 



+ 



P 



P 
P 



7 7* 



dZ 

d"fin 

dZ 

n t: 

O-fin 

dZ 



(2.22) 
(2.23) 
(2.24) 



where the notation f and P stand for cut and uncut Pomeron respectively and only second order terms in 7 ("fin) 
are shown. The next step it to use AGK cutting rules to write the resulting evolution equation 



dZ_ 
dY 



r(i 



37 37 
2 ) (7 - 7 2 ) 7r + T(l - 2)(2 7 2 - 4 77m + jf n + lin ) — 
07 07™ 



(2.25) 



It can be easily seen that the second term reproduces the first term for 2 7 = 7 i„ . We note that u = 1 — 7 and 
v = 1 — jin and thus Eq. ( 2.25| ) for u and v reads as 



«_-r(i~ a Mi-)f£-r(i 



dZ 

2) (2m 2 - Auv + v 2 + v) — 
ov 



(2.26) 



The description in terms of generating function becomes clear for Markov chain. We use the definition of 
generating function given by Eq. (2.14) to find the differential equation for probabilities 

dP (n,m\Y) 



d(T(l -» 2) Y) 
(P —> P + P) 

(f -> r + f ) 
(f -> r + p) 
if -> p + p) 



- nP(n,m\Y) + (n - 1) P (n - 1, m\Y) 
+ mP(n,m\Y) — (m — 1) P (n, m — 1\Y) 

- 4m P (n,m\Y) + 4mP (n - l,m\Y) 

+ 2 m P(n,m\Y) - 2 (m + 1) P (n - 2, m + 1|V) 



(2.27) 
(2.28) 
(2.29) 
(2.30) 



where each line describes specific process of Pomeron splitting discussed above. It is instructive to compare this 
equation to Eq. (2.6). Each of Eq. ( 2.27 ) - Eq. ( 2.30 ) consists of two terms: one describes to increase of probability 
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to find n-Pomerons due to decay of one Pomeron to two and one is responsible for the decrease of this probability 
since one of n Pomerons can decay. In all equations, except Eq. ( |2.28|) and Eq. ( 2.30| ), the increase leads to the plus 
sign and a decrease to the minus sign. However, in Eq. ( 2.2S ) and Eq. (|2.30|) the signs are opposite in accordance 
with the AGK cutting rules (see Fig. 0). In this case we have to say that decay P — ► P + P and P — > P + P 



have negative amplitudes. It should be stressed that in terms of the amplitude (see Eq. (2.25)) we obtain positive 
cross sections with different multiplicity of produced particles. 



Equations Eq. ( 2.27 ) -Eq. ( 2.30 ) give clear probabilistic interpretation of all Pomeron splitting processes under 
discussion. 



It turns out that Eq. ( 2.26| ), being a typical Liouville equation, has the solutions that depends on two variables 



6 = T(l -> 2) Y + In and & = T(l -» 2) Y + In ■ 



The general solution of this equation is given by an 



l-(2u-v) ■ 

arbitrary function of variables £i and £2, which in our case can be written as sum of two functions, namely, 



Z = F 1 iT(l^2) Y + ln 



1-u 



F 2 lT(l^2) y + in 



2u 



1 - (2u-v) 



For our initial condition Eq. (2.20) the solution reads 

2ue -n^2)Y 



z = 



(2u — v)e 



-r(i^2)y 



1 + u{e- T ^ 2 ) Y - 1) 1 + (2u - i;)(e- r ( 1 - 2 ) Y - 1) 



(2.31) 



(2.32) 



One can easily see that this solution satisfies both initial Eq. ( 2.20| ) and boundary Eq. ( 2.1 5| ) conditions 



2.4 The generating functional for the multiparticle production: non-linear equation 



Using our initial condition Eq. (2.20) we can rewrite linear differential equation Eq. (2.26|) as non-linear one. We use 



the fact mentioned above, namely, that Eq. (2.26) is differential equation of two variable £1 and £2 and, thus, has no 
separate dependence on Y. This means that the differential equation Eq. ( 2.26Q written at some rapidity Y keeps 
the same form for any rapidity. We pick initial rapidity Y — and substitute generating function given by Eq. (2.20) 



into Eq. fl2.26|) 



dZ_ 
dY 



- r(l -> 2)(2u 2 - Auv + v 2 + v) 



(2.33) 



Now use initial condition from Eq. fl2.g|) for generating function with no cut Pomerons. In terms of Zq and Z 



Eq. Q2.33D reads 

dZ 

— = -T(1^2)(2Z 2 -4Z Z + Z 2 + Z) (2.34) 
We identify scattering amplitude and diffractive cross section with iVo(7|Y) = 1 — Z (l — j\Y) and N(^,^i n \Y) = 



1 - Z(l -7,1 - 7 m |F), respectively. For N (-y\Y) and N(~f,~/ in \Y) Eq. (|2.34|) reads 
dN{ lnm \Y) 



dY 



= -r(l - 2) {27V 2 ( 7 |r) - 4iV (7in^(7,7m|y) + N 2 ( 1 , lvn \Y) + N{ 1 , lln \Y)} 



(2.35) 



This equation has the same form as the equation for the diffraction production 1 obtained by Kovchegov and Levin 
(see Eq. ( |2.3| )), but it is written for a general functional and describes not only diffractive production but also the 
processes of particle production with any value of multiplicity. The fact that Eq. ( |2.3| ) has a simple generalization 
on the general case of QCD gives us a hope to generalize this equation. 



lr The difference in overall minus sign corresponds to different definitions of rapidity variable moving in opposite direction. 
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2.5 The generating functional for the multiparticle production: consistency with the AGK cutting 
rules 



We want to check the consistency of our solution Eq. (2.32) with the AGK cutting rules in an explicit way. To do 
this we calculate cross section of a process with a given multiplicity k from both generating function given by 
Eq. (2.32) and directly from the AGK cutting rules, and compare them. We define cross section with multiplicity k 

as 



(fc) _ 1 d*N(~,, 7in \Y) 

u — - • ~ ' hin=0 Jin 



(2.36) 



k\ 87* 

where AT( 7 , 7in |F) = 1 - Z{1 - 7 , 1 - 7in |F). 

As an example we pick multiplicity to be that of one cut Pomeron, i.e. k = 1 for any number of uncut Pomerons. 
From Eq. (2.36) with k = 1 we get 



line 



r(i^2)r 



(l + 2 7 (e r ( 1 ^ 2 ) y -l)) 2 

On the other hand we can use coefficients for multiple Pomeron exchange from the AGK cutting rules 
case cross section for multiplicity of k cut Pomerons reads 

oo 

a« = ^(-l)™- fe C£(2 7 )"e™ r ( 1 ^ 2 ) Y 

n=Q 

where e 1 ^ 1 "* 2 ^ stands for Pomeron propagator. 



(2.37) 
In this 

(2.38) 



For k — 1 we sum over n in Eq. ( 2.38 ) and obtain 



r (i) 



2 7 e 



r(i^2)F 



(1 + 2 7 e r ( 1 ~> 2 ) 5 2 



(2.39) 



In high energy limit 7 ; n = 2 7 and e r ^ 1 ^ 2 ^ Y — 1 ~ e r ^^ 2 ^ Y , thus Eq. (2.37) reproduces Eq. ( [2.39 ). It can be 
easily shown that this holds for any value of k. 

We can further compare cross sections obtained from generating function and direct summation of fan diagrams 
using AGK rules |l0| . In our approach this corresponds to 

a 8d = N(rf, j m = 0\Y) a m = A( 7 , 7m = 2 7m |T) (2.40) 

where iV( 7 , 7 i„|Y~) = 1 — Z(l — 7 , 1 — 7 m|T). The resulting expressions are identical and given by 



2 7 2 e r(i->2)i- (e r(i^2)Y _ 1} 



(1 + 7 ( e r(i^2)y _ + 27 ( e r(i^2)y _ ^ 



2 7 e 



r(i-)-2)r 



l + 2 7 (e r ( 1 - 2 ) y -l) 



(2.41) 



2 7 e 



r(i^2)Y 



1 + 7 (eT(i-f2)y _ x) 
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3. Pomeron loops 



3.1 Evolution equation with loops 



Now we want to account for contributions of Pomeron loops. This problem has already been solved for a case with 
no cut Pomerons |36) . The master equation for no cut Pomerons is given by 



«— ra^nwi-og + ^i>(. (3-D 



Unfortunately, one cannot generalize Eq. (2.26) for diffractive processes by just adding second order derivative terms 
by analogy with Eq. (3.1). The reason for that is because this type of equation would include diagrams that does 
not exist. For example, diagrams of the type f — > P + P — > P where cut Pomeron splits to two uncut Pomerons 
with further merging to uncut Pomeron are forbidden. To resolve this problem we introduce two separate variables 
w and w for uncut Pomeron being in amplitude or conjugate amplitude, respectively. Naturally, these two subsets 
of Pomerons evolve separately till the cut Pomeron is introduced, alternatively, their evolutions mix only via cut 



Pomeron. Using our previous discussions we can readily write this new type of evolution equation based on Eq. (2.26) 



dZ . . f . ,dZ , s dZ 

dZ 

- T(l -> 2)(2ww - 2wv - 2wv + v 2 + v)— (3.2) 

ov 

+ »rfl~i){o-^)g + 

- ¥ {2 1] { 2{v - "*»^ + 2( * - + 2( " " ^iSr, + <" ^ " 2) S > z 



where the first and second lines are generalization of Eq. ( |2.2q ) for variables w and w, the third line corresponds to 
Eq. ( |3.l|) and the last line a little bit more explanation. We have to introduce a new generating function 



Z(w, w, v\Y) = ^ X] zZ P ( fc ' l > ™\Y)w k w l v m (3.3) 

k=0 1=0 m=0 



where P(k,l,m\Y) stands for probability to find k uncut Pomerons in the amplitude, I uncut Pomerons in the 
conjugate amplitude and m cut Pomerons at some rapidity Y . For the last term of Eq. ( |3.2j ) we write Markov chain 
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in a similar manner we did it for Eq. (2.27)-Eq. (2.3C), namely, 



dP(k,l,m 


1*0 








dY 








(P - 


-> P - 


h P) 


+ r(i -> 


2) {(fc - l)P(fc - 1, I, m\Y) - kP(k, I, m\Y)} 


\ir - 


f - 


- P) 


+ r(i 


2) {(I - l)P(fc, I - 1, m|F) - ZP(fc, Z, m|y)} 


(? 


- f 


f f) 


- r(i- 


2) {(to - l)P(fc, Z, to - l\Y) - mP(k, I, m\Y)} 


(f - 


-> P - 


h f ) 


+ 2T(1- 


-> 2) {mP(fc - 1, /, m|Y) - mP(fc, I, m\Y)} 


(r - 


-> P - 


" f ) 


+ 2T(1- 


-> 2) {wP(fc, Z - 1, to|F) - mP(fc, Z, m|F)} 


(r 


-> P 


f P) 


- 2T(1- 


-> 2) {(to + l)P(fc - 1, 1 - l,m + l|y) - mP(fe, Z, m|y)} 


(P 


■f P- 


+ P) 


+ ^r(2 


— > l) {fc(fc + i)p(fc + 1, z,to y) - fc(fc - i)p(fc, z,m|y)} 


(p - 


f P- 


> P) 


+ \?<? 


-> 1) {Z(Z + l)P(fc, l + 1, m|y) - 1(1 - l)P(fc, z, m|y)} 


(f 


ff- 


♦ P ) 


- 2^ 


-> i) {to(to + i)p(fc, z, to + i|y) - to(to - i)P(k, z, m|y)} 


(f 


■f P- 


♦ P) 


2 

- i r < 2 


-> 1) {(fc + l)mP(fe + 1, Z, m|y) - fcTOP(fc, Z, m\Y)} 


(r - 


f P- 


> f ) 


- 5 r < 2 


-> 1) {(Z + l)mP(k, I + 1, m|y) - ZTOP(fc, Z, m|y)} 


(p - 


f P- 


> f ) 


2 

- - r(2 

2 V 


-> 1) {(fc + 1)(Z + l)P(fc + 1, Z + 1, m\Y) - fcZP(fc, Z, m\Y)} 



(3.4) 



where P denotes Pomeron in conjugate amplitude and factor of ^ accounts for a fact that the Pomerons are identical 
in this approach. Each line in Eq. (3.4) as in Eq. (2.27)-Eq. (2.30) has a clear probabilistic interpretation. We 
multiply Eq. (EOI) by w k w L v m , sum over all fc, I and to, and using definition Eq. (|3.3|) obtain Eq. (|3.2|). 



Even in our simple model with no coordinate dependence Eq. ( p.2| ) is too much complicated and its solution is 
still to be found. But before the solution is found we can see that putting w — w = u, i.e. making no difference 
between uncut Pomerons, first two lines of Eq. ( p.2j ) correctly reproduce Eq. ( 2.26| ). Moreover, we can check it further 
and using initial condition Z(w, w, v\Y = 0) = v we can perform iterations 



= -r(i 

dY v 



2)(2i 



2wv — 2wv + v 2 + v ) 



giving 



Zi 



-r(l -> 2)(2wZi - 2tro - 2ii)iJ + w 2 + u)y 



(3.5) 



(3.6) 



At the next step 



y 2 

T 2 (l -> 2) {tu(l - w)(2u) - 2u) - w(l - u))(2w - 2u)} — 

y 2 

T 2 (l -> 2)(2uw - 2wv ~ 2wv + v 2 + v)(-2(w + w) + 2v + 1)— (3.7) 

l y 2 

-r(l -> 2)r(2 -> 1) {2(w - wiD)(-2) + 2(v - ww)2 + 2(« - wv)2 + (v - v 2 )(~2)} — 



It is easy to see from Eq. (3.7) that Eq. ( |3.2|) correctly reproduces the sign and combinatorics coefficient of Pomeron 
loops in accordance with AGK cutting rules. For example, in the first term in the third line — — * 2)T(2 — > 
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Figure 5: The renormalization procedure in the case of the simplest 'fan' diagram. 



l)2(i> — ww)(— 2), the term proportional to v describes the Pomeron loop of the type f — > P + P — > ^ . This loop 
has factor of 2 and brings plus sign as expected. Similarly, the terms proportional to v in the second and third term 
of the last line correspond to loops of the type f — ► f + P — ► f and f —> f + P —> respectively. Each of 
them brings factor of —2, putting P = P we have a factor of —4 which is in a agreement with AGK cutting rules. 



Thus, we expect Eq. (3.2) properly include Pomeron loops into evolution. 

At the end of this part we would like to mention that the problem has been already treated by Ciafaloni and 
Marchesini EJ many years ago, but in terms of RFT Lagrangian. They separated uncut Pomerons in amplitude 
and conjugate amplitude by introducing different fields (0+ and <f>-). It is interesting to notice that, just like in 
our case, the authors also have partial diagonalization of Lagrangian in terms of new variable <j) + + 0_ — i<f> c . By 
partial diagonalization one should understand diagonal vertices only in one rapidity direction. In our case this can 
be formulated as separation of evolutions of Pomerons corresponding to w, w and w + w — v in MFA. 



3.2 A new method of summation of the Pomeron loops: improved Mueller-Patel-Salam-Iancu ap- 
proach 



Solution to Eq. (3^) will give the generating functional. However, Eq. (3^2 ) is based on the probabilistic interpretation 
of the BFKL Pomeron calculus in terms of Markov process. In this subsection we would like to suggest a different 
interpretation with a different procedure of summation of the Pomeron diagrams. 

Firstly we consider the simplest 'fan' diagram of Fig. |^. It can be calculated in an obvious way, namely, 

A (Fig- H) = 7G(Y-0) - A 7 2 / dy 1 G(Y-y 1 )G 2 (y 1 -Q) (3.8) 

Jo 

= 7 e Ay - A 7 2 C dy x e A ( y +^) (3.9) 



7 e 



A 6 



„A Y 



7 2 e 2Ar + ( 7 + 7 2 )e Ay 



7 e + 7«e 



where T(l — > 2) (see Eq. (2.7)) is denoted as A , 7 is the amplitude of the Pomeron interaction with the target and 
G(Y — y) stands for the Green function of the Pomeron G(Y — y) = exp (A (Y — y)). 

As one can see, the integration over yi reduces the diagram in Fig. ^| to two contributions: the exchange of two 
non- interacting Pomerons and the exchange of one Pomeron with the renormalized vertex: 7^ = 7 + 7 2 . In Fig. || 
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Figure 6: The renormalization procedure in the case of the 'fan' diagram of the second order. 



is shown the Pomeron 'fan' diagram of the second order which we have to integrate over two rapidities y\ and y 2 . 
The result is 



A(Fig.§ = 
= 2A 2 7 3 
= 2A 2 7 3 



dyi 



dyi 



Vl dy 2 G(Y ~ Vl ) G( Vl - 0) G(y 2 - 0) G 2 (y 2 - 0) 
v 1 ... / 1 



(3.10) 



dy 2 e 



A(Y+ Vl +y 2 ) 



2 „,3 



7' 







3Ar 



2 AY 



AY 



2A 2 A 2 

Adding the contributions of this diagram and the diagrams of Fig. || we obtain 
A {Fig. § + A{Fig.% = 



2 A 2 



7 3 e 3AY - 7(7 + 7 2 ) e 2AY + (7 + 7 2 + 7 3 )e 3Ar 



7 3 e 3AY 



277^ e 2A - 



1RG 



(3.11) 

ay 



Therefore, one can see that the scattering amplitude can be rewritten as exchanges of the Pomerons without in- 
teraction between them but with Pomeron - particle vertex. In the dipole model this vertex has a meaning of the 
amplitude of two dipole interaction in the Born approximation of perturbative QCD. 

These two examples illustrates our main idea: the BFKL Pomeron calculus in zero transverse dimensions can be 
viewed as the theory of free, non-interacting Pomerons whose interaction with the target has to be renormalized. It 
is easy to see that in the MFA we can rewrite the master equation (see Eq. (2.7) and Eq. (2.11)) in the form 



dN«( lR \Y) dN«{ lR \Y) 
= 1 (1 — ► I) 7^ 



dY 



9jr 



with 



1R = 



(3.12) 



(3.13) 



The way Eq. ( 3.13] ) has started to build in perturbation expansion we have shown in Eq. (3.S) and Eq. (3.10). 

The general solution of Eq. (3.12) is the system of non-interacting Pomerons and the scattering amplitude can 
be found in the form 



No(lR.\Y) = ]T C nl n R G n {Y-Q) 



(3.14) 



n=l 
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where coefficients C n could be found from the initial conditions, namely, from the expression for the low energy 
amplitude. In particular, the initial condition 



iVo(7fl|^ = 0) = 7 = 7fl/(l+7fl) 
generates C„ = (—1)™ and the solution is 



N (~/r\Y) = 



AY 



1 



AY 



(3.15) 



(3.16) 



The initial condition of Eq. ( |3.15 ) has a very simple physics behind it. It describes the independent (non-correlated) 
production of the Pomerons at low energy but with the only condition that one Pomeron lives shorter time than the 
second one (see Fig. |?]). If all n Pomerons were emitted by the same dipole (see Fig, ^ -b) this condition leads to 
Glauber factor 1/n! leading to Nq^rIY = 0) = 7 = 1 — exp(— 7^). instead of Eq. ( 3.15 ). However, if Pomerons are 
produced as the consequent decays (see Fig. Q-a) the factor is equal to 1. (—1)" comes from the Glauber screening, 
resulting in Eq. fl3.15p . In QCD we have strong evidence that the second case is correct [M . 

The analysis of enhanced diagrams we start from the first dia- 

Irl ^ > * 

gram of Fig. H. It leads to the following contribution 



A (Fig. § 
-A 2 7 2 



dyi 



dy 2 G(Y -y 1 )G 2 (y 1 -y 2 )G(y 2 



(3.17) 

-o) 



i 2 7 2 



dyi 




dy 2 e 



A(Y+ V i-y 2 ) 



c3 



g 




3 0608 
3 


5606® 

O 
O 
O 






O 

P 






P 





2 2 AY 1 2 

7 e + 7 e 





AY 



A Ye 



AY 



Figure 7: The diagram that illustrates the initial 
condition of Eq. (3.15). 



where T(2 ->■ 1) = A7 2 (see M). 



Adding Eq. (3.17) to the exchange of the one Pomeron we obtain that the exchange of one Pomeron and the 
enhanced diagram of Fig. @ can be written in close form 



One Pomeron exchange + A (Fig. ||) = 7^ e 



An Y 



7 2 e 2Ay 



with 



,(2) 



(3.18) 



(3.19) 



7-R = T"' = 7 + 7; A fl = A + 7A; 

It is easy to see the Eq. (|3.17 ) can be viewed as the expansion to first order of Eq. ( 3.1 8| ) 

Therefore, the Pomeron loops can be either large (of the order of Y) and they can be considered as un-enhanced 
diagrams ; or small (of the order of 1/A) and they can be treated as the renormalization of the Pomeron intercept. 

In QCD A oc as while 7 oc a 2 s . Therefore, the renormalization of the Pomeron intercept A is proportional 
to a s . We can neglect this contribution since (i) there a lot of a 2 s corrections to the kernel of the BFKL equation 
that are much larger than this contribution; and (ii) in the region of Y ^ l/a 2 s 1 where we can trust our Pomeron 
calculus (see introduction) (Ar — A) Y -C 1. 

Concluding this analysis we can claim that the BFKL Pomeron calculus in zero dimensions is a theory of non- 
interacting Pomerons with renormalizcd vertices of Pomeron-particle interaction. In the dipole language, it means 
that we have a system non-interacting Pomerons with a specific hypothesis on the amplitude of the dipole interactions 
at low energy. For the problem that we are solving here, namely, when we have one bare Pomeron at low energy, 
this amplitude is determined by Eq. ( 3.1 5| ) . 
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For such a system we can calculate the scattering amplitude using a method suggested by Mueller, Patel, Salam 
and Iancu and developed in a number of papers (see |33], 0, 0, |39], [|, [| and references therein). This method 
suggests that the scattering amplitude can be calculated using the i-channel unitarity constraints which is written 
in the following way (assuming that amplitudes at high energy are pure imaginary, N — Im A) : 

N([...]\Y) = N([...]\Y-Y';P^nP)(g)N([...]\Y';P^nP) (3.20) 

where (^) stands for all needed integrations while [. . . ] describes all quantum numbers (dipole sizes and so on ) . 
The correct implementation of this leads in our case to the following formula (see also |39| [l2| ^| ) 

N MPS I{Y) = ! _ exp |_ 7 B^^^|^MFA (7i|y _ FOiv M^ (7 2 |FOU = 7| = o (321) 

where J\[ MFA (Y, jr) is given by Eq. ( |3.16j )(see also Eq. ( 3.12 )) in the mean field approximation and j BA oc a 2 s is 
the scattering amplitude at low energies which is described by the Born approximation in perturbative QCD. The 
difference of Eq. (3.21) from the original MPSI approach is the fact that this equation does not depend on the value 
of Y' and, because of this, we do not need to choose Y' = Y/2 for the best accuracy. 



Substituting Eq. (3.16) in Eq. (3.21) we obtain 

pjMPSi ( 7 BA {Y) = l cxp ^ BA J {1 ^ )Y T ( 0, „, BA J (1 ^ 2)Y ) (3.22) 

T (0, x) is the incomplete gamma function (see formulae 8.350 - 8.359 in p3[). 



We claim that Eq. ( 3.22 ) is the solution to our problem. One can easily see that Nq ("f\Y) — > 1 at high energies 
in contrast to the exact solution with Hamiltonian of Eq. ([O]). The exact solution leads to the amplitude that 
vanishes at high energy (see [Q, ||). As have been mentioned the solution depends crucially on the initial condition 
for the scattering amplitude at low energies. For Eq. (3.22) this amplitude is equal to 

oo 

N$ IPSI ( 7 |Y - 0) - (-^ n+1 n[ b BA T ( 3 - 23 ) 

71 = 1 

with j BA oc a 2 s . This equation reminds us the ultraviolet renormalons contribution and calls for better understanding 
of the non-perturbative contribution. 



We can rewrite Eq. (3.21) in more convenient form using the Cauchy formula for the derivatives, namely, 



d n Z MFA( lR \ Y ) .If Z MFA {i R \Y) 



■dj' R ; (3.24) 



Contour C in Eq. (3.24) is a circle with a small radius around 7^ = 0. However, since function Z does not grow at 
large 7^ for n < 1 we can close our contour C on the singularities of function Z . We will call this new contour Cr. 

n mps I{y) = 1 exp{- 7 S ^^} N^ A ({ lR \Y-Y>) (( 1r \Y>) | 7 > =7 * =0 

v BA\ n 1 r yMFAu.l W _ V>\ r 7 MFA^,2 



7_) 1 I ^ Z™r\ lR \Y-Y>) I ^ Z m ^( 7r \Y') 

n\ 



1 - V V ' > n\n\ — & d lR - //f' ' J> d lR ., .. 



71=1 

' / / d%dj R f / 7^7| \ IrIr t ( q Inll 



(2TTi) 2 J J 7^ 7| \ \^^BAgT(l >2)Y J ryBAgT(l >2)Y ^ ' ryBA^T(l >2)Y 

Z MFA ( -l ) Z MFA ( -2) ( 3 .25) 
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Figure 8: The renormalization procedure in the case of the simplest enhanced diagram. 

Here we introduce new variables 7^ = 7^ exp (T(l — > 2) (Y — Y')) and j R = j R exp (r(l — > 2) Y - '). In these new 
variables 

Z MFA (?r) = Y^l-T ; gMFA m = _L_ (3 26) 



Closing the integration on the poles -yjj = — 1 and 7^ = —1 we obtain the formula of Eq. ( |3.22 ) 



3.3 The generating functional for the multiparticle production with Pomeron loops 

Using the result of the previous section we will derive the formula in the MPSI approach for the general functional, 



defined by Eq. ( p.3[) . This formula is based on the solution of Eq. ( |3.2|) but without the secondary derivatives. Such 
a solution gives the MFA approximation to our problem and we denote it as Z MFA (w,w,v\Y). The equation for 
Z MFA (w, w, v\Y) looks as follows 

Z(w,w,v\Y) = (3.27) 
we -r(i^2)Y ^ e -r(i^ 2) y ( w + ^_ i; ) e -r(i^2)y 

1 + w{e- r ^-' 2 ) Y - 1) + 1 + w{e-n^)Y _ ~ l + ( w + w - w)( e - r ( 1 ^ 2 ) Y - 1) 



Using the renormalized 7 of Eq. (p,13j) we can rewrite Eq. ( 3.27 ) in a different form, namely. 



Z MFA { lR ^ R , lin , R \Y) = ' ' ' 



1 + 1R e r ( 1 ^ Y l + jne^ 1 ^ 2 ^ l + ^e r ( 1 ^ 2 ) y 
1 1 1 



1 + 1R 1 + 1R 1 + £,R 

where we use a notation £ = 1 — w — w + v = 7 + 7 — 7™ and 



(3.28) 



Zr = ir + ir - n n ,R = r^j ; f = ; (3.29) 
j R = lR e rd-2)y . = fl = ^ e r ( i^ 2) y . ^ = ^ e r(i^2) y . 



The first of Eq. ( 3.29| ) is the definition of ji n ,R- 
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The general formula for the amplitude in the MPSI approach has the form 



N 



MPSI ^BA^BA^ 



exp 



tMFA 



,,BA 



d d 



RA 



d d 



RA 



d 



d 



dlh.dl 2 R ' dj R dj R + lin d^ R d^ R 

tMFA 



(3.30) 



{iWr,iL,r\Y ~ Y') Z MFA {l R ,f R nl,R\Y') | 7i = ^1=^^ ^ o 



exp 



■7 



BA 



d d 



7 



BA 



d d 



+ 7 



BA 



d 



d 



V L 9% ff? R ' d% d~f R ' '«* 9t4,* d^ n>R 

7 MFA (~1 =1 -1 \ ryMFA ( ~2 = 2 -2 \ I 

z l7ji,7ii)7m,iij z (,7fl I 7ij J 7 4 ™,flj 1^=73=71=73=71 =72 = a 



1 



where 7 s " 4 and 7^ are the elastic and inelastic amplitudes of interaction of two dipoles at low energy which are 
calculated in QCD in the Born approximation. The plus sign in Eq. (|3.30D in front of jg A reflects the 

fact that the Pomeron loop with two cut Pomerons does not have a negative contribution unlike in the case of uncut 



Pomerons. The sign between the exponent and unity in Eq. (3.30) could be easily checked noticing the first term of 
the expansion of the exponent correctly reproduces the positive contribution of ji n ( cut Pomeron). 

The nice feature of this equation that one can see that the result does not depend on the value of an arbitrary 
chosen rapidity Y'. Using the explicit form for Z MFA given by Eq. ( 3.28 ) we can calculate Z MPSI in closed form 
denoting by 

G(x) = exp(-) -r(0, -) (3.31) 



N 



MPSI 



is equal to 



N 



MPSI 



(7 



BA 



1 BA 

1 in 



Y) 



(i- G ( 



7 _ba e r(i^2)y 



1 - G (2 7 



„BA 



7f/) 



T(i-2)y 



(3.32) 



The useful formulae for getting Eq. (3.32) are the following 

Qk gh Ql 2 1 



(-1)' 1+ ' 2 (h +h + k)\ 



dlin d ~f h ®l l2 1 + 7 + 7 - lin 

and equations 8.350 - 8.359 for the incomplete gamma function T (0,x) in 



(3.33) 



Eq. ( 3.32 ) allows us to calculate the cross section with fixed multiplicity of produced particles. Namely, the cross 
section for the processes with k < n > particles in the final state, where < n > is the mean multiplicity in our 
reaction, can be calculated as 



<Tk(X) = 



1 



Qh: 



■ W/Wl U,B 1 , B . 1 



7 ? 7i'i 



\Y) 



(7f/ = 2 7 SA ) 



(3.34) 



Here we use that 7 



BA 



27 in the Born approximation of QCD. 



It is interesting to check the general equation (see Eq. (3.32)) calculating two known cases: the diffractive 



dissociation process and the total inelastic cross section. The first one can be calculated using Eq. (3.32) with 
7^? = 0. The answer is 



x^Iff 1 {i ba \y) = 2iV MPS7 (7 SA |r) - 7v MPS/ (2 7 SA |r; 



(3.35) 



where N is given by Eq. ( 3.22 ). Eq. ( 3.35 ) is a direct consequence of the unitarity constraints( see Eq. ( 2.16 )). As 
you can see from Eq. (2.32) the same formula determines the diffractive production in the mean field approximation. 
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The value of the total inelastic cross section, that is equal the sum of diffractive production and inelastic cross 
section, stems from Eq. ( |3.32| ) for 7^* = 2 j BA and it is equal to 

N^FLlasHc h BA \y) = 2N MPSI (l BA \Y) = Nfff ( 7 *» + N& ai (l BA \Y) (3.36) 

which is actually the unitarity constraint itself ( see Eq. (|2.16| )). One can see that the inelastic cross section is 
determined by 

N^ SI (l BA \Y) = 1 - G( 2l BA e^ Y ) (3.37) 



4. Conclusions 



In this paper we introduce a new generating function for the processes of multiparticle production both for the mean 
field approximation (see Eq. ( 2.14 ) and for the general case (see Eq. (3.3)). For general case where the Pomeron loops 
has been taken into account, we obtain the linear evolution equation for the generating function (see Eq. (3^2)) while 
in the mean field approximation we proved both the linear evolution equation (see Eq. ( 2.26| )) and the non- linear 
equation (see Eq. (2.34)). The last one is the generalization of Kovchegov-Levin equation for diffractive production 
to a general case of the processes with arbitrary multiplicities. Since this equation is proven for the general QCD 
case we hope that the equations for the general generating function can be proven for the real QCD evolution. 

The second result of the paper is the new method of summing the Pomeron loops. We argued that the sum 
of all Pomeron diagrams, including loops, in the kinematic region of Eq. (1.2) can be reduced to the diagrams 
of the Pomeron exchanges without interactions between Pomeron if we renormalize the amplitude of low energy 
interaction. Based on this result we suggest a generalization of the Mueller-Patel-Salam-Iancu method of summation 
of the Pomeron loops. In particular, we calculated the new generating function for the inelastic processes in the 
improved MPSI approximation ( see Eq. ( 3.32j )). 

We would like to stress that we firmly believe that the scattering amplitude, calculated using this method, leads 
to a correct answer to the old problem: the high energy asymptotic behaviour of the scattering amplitude at ultra 
high energies beyond of the BFKL Pomeron calculus in kinematic region of Eq. (1.2). 

We hope that both results will lead to new simplifications in the case the BFKL Pomeron calculus in QCD (in 
two transverse dimensions) . The general case of the BFKL Pomeron calculus in QCD will be addressed in a separate 
paper. We would like also to mention that this case has been started to discuss in p4j. 



Acknowledgments: 



We are very grateful to Jochen Bartels, Asher Gotsman, Edmond Iancu, Michael Kozlov, Uri Maor and Al Mueller 
for fruitful discussions on the subject. Special thanks of A. P. goes to Jochen Bartels for his hospitality and support 
during his visit to DESY. This research was supported in part by the Israel Science Foundation, founded by the 
Israeli Academy of Science and Humanities and by BSF grant # 20004019. 

We thank the Galileo Galilei Institute for Theoretical Physics for the hospitality and the INFN for partial support 
during the completion of this work. 



References 

[1] D. Amati, L. Caneschi and R. Jengo, Nucl. Phys. B 101 (1975) 397 V. Alessandrini, D. Amati and R. Jengo, Nucl. 
Phys. B 108 (1976) 425; R. Jengo, Nucl. Phys. B 108 (1976) 447; D. Amati, M. Le Bellac, G. Marchesini and 



- 18 - 



M.Ciafaloni, Nucl. Phys. B112 (1976) 107; D. Amati, G. Marchesini, M.Ciafaloni and G. Parisi, Nucl. Phys. B114 
(1976) 483; M. Ciafaloni, M. Le Bellac and G. G Rossi, Nucl. Phys. B 130 (1977) 388. 

[2] E. Levin, Nucl. Phys. A763, 140 (2005), arXiv:hep-ph/0502243. 

[3] K. G. Boreskov, A. B. Kaidalov, V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C 44 (2005) 523 
[arXiv:hep-ph/0506211]. 

[4] A. I. Shoshi and B. W. Xiao, Phys. Rev. D 73, 094014 (2006) [arXiv:hep-ph/0512206]. 

[5] M. Kozlov and E. Levin, Nucl. Phys. A 779 (2006) 142 [arXiv:hep-ph/0604039]. 

[6] A. I. Shoshi and B. W. Xiao, arXiv:hep-ph/0605282. 

[7] J. P. Blaizot, E. Iancu and D. N. Triantafyllopoulos, arXiv:hep-ph/0606253. 

[8] S. Bondarenko, L. Motyka, A. H. Mueller, A. I. Shoshi and B. W. Xiao, arXiv:hep-ph/0609213. 

[9] M. Kozlov, V. Khachatryan,E. Levin and J. Miller, arXiv:hep-ph/0610084. 
[10] L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rep. 100, 1 (1983). 
[11] A. H. Mueller and J. Qiu, Nucl. Phys. ,427 B 268 (1986) . 

[12] L. McLerran and R. Venugopalan, Phys. Rev. D 49,2233, 3352 (1994); D 50,2225 (1994); D 53,458 (1996); D 59,09400 

(1999) . 

[13] I. Balitsky, [arXiv:hep-ph/9509348]; Phys. Rev. D60, 014020 (1999) [arXiv:hep-ph/9812311] Y. V. Kovchegov, Phys. 
Rev. D60, 034008 (1999), [arXiv:hep-ph/9901281]. 

[14] J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert, Phys. Rev. D59, 014014 (1999), [arXiv:hep-ph/9706377]; 
Nucl. Phys. B504, 415 (1997), [arXiv:hep-ph/9701284]; J. Jalilian-Marian, A. Kovner and H. Weigert, Phys. Rev. 
D59, 014015 (1999), [arXiv:hep-ph/9709432]; A. Kovner, J. G. Milhano and H. Weigert, Phys. Rev. D62, 114005 

(2000) , [arXiv:hep-ph/0004014] ; E. Iancu, A. Leonidov and L. D. McLerran, Phys. Lett. B510, 133 (2001); 
[arXiv:hep-ph/0102009]; Nucl. Phys. A692, 583 (2001), [arXiv:hep-ph/0011241]; E. Ferreiro, E. Iancu, A. Leonidov 
and L. McLerran, Nucl. Phys. A703, 489 (2002), [arXiv:hep-ph/0109115]; H. Weigert, Nucl. Phys. A703, 823 (2002), 
[arXiv:hep-ph/0004044] . 

[15] E. Levin and K. Tuchin, Nucl. Phys. A693 (2001) 787 [arXiv:hep-ph/0101275]; A691 (2001) 779 
[arXiv:hep-ph/0012167]; B573 (2000) 833 [arXiv:hep-ph/9908317]. 

[16] E. Iancu, K. Itakura and L. McLerran, Nucl. Phys. A708 (2002) 327 [arXiv:hep-ph/0203137]. 

[17] N. Armesto and M. A. Braun, Eur. Phys. J. C20, 517 (2001) [arXiv:hep-ph/0104038] ; M. Lublinsky, Eur. Phys. J. 
C21, 513 (2001) [arXiv:hep-ph/0106112]; E. Levin and M. Lublinsky, Nucl. Phys. A712, 95 (2002) 
[arXiv:hep-ph/0207374]; Nucl. Phys. A712, 95 (2002) [arXiv:hep-ph/0207374]; Eur. Phys. J. C22, 647 (2002) 
[arXiv:hep-ph/0108239]; M. Lublinsky, E. Gotsman, E. Levin and U. Maor, Nucl. Phys. A696, 851 (2001) 
[arXiv:hep-ph/0102321]; Eur. Phys. J. C27, 411 (2003) [arXiv:hep-ph/0209074]; K. Golec-Biernat, L. Motyka and 
A.Stasto, Phys. Rev. D65, 074037 (2002) [arXiv:hep-ph/01 10325]; E. Iancu, K. Itakura and S. Munier, Phys. Lett. 
B590 (2004) 199 [arXiv:hep-ph/0310338]. K. Rummukainen and H. Weigert, Nucl. Phys. A739, 183 (2004) 
[arXiv:hep-ph/0309306]; K. Golec-Biernat and A. M. Stasto, Nucl. Phys. B668, 345 (2003) 

[arXiv:hep-ph/0306279]; E. Gotsman, M. Kozlov, E. Levin, U. Maor and E. Naftali, Nucl. Phys. A742, 55 (2004) 
[arXiv:hep-ph/0401021]; K. Kutak and A. M. Stasto, Eur. Phys. J. C41, 343 (2005) 
[arXiv:hcp-ph/0408117]; G. Chachamis, M. Lublinsky and A. Sabio Vera, Nucl. Phys. A748, 649 (2005) 
[arXiv:hep-ph/0408333]; J. L. Albacete, N. Armesto, J. G. Milhano, C. A. Salgado and U. A. Wiedemann, Phys. Rev. 
D71, 014003 (2005) [arXiv:hep-ph/0408216]; E. Gotsman, E. Levin, U. Maor and E. Naftali, Nucl. Phys. A750 (2005) 
391 [arXiv:hep-ph/0411242]. 

[18] A. H. Mueller, A. I. Shoshi and S. M. H. Wong, Nucl. Phys. B 715, 440 (2005) [arXiv:hep-ph/0501088]. 

[19] E. Levin and M. Lublinsky, Nucl. Phys. A 763, 172 (2005) [arXiv:hep-ph/0501173]. 



- 19 - 



[20] E. Iancu and D. N. Triantafyllopoulos, Nucl. Phys. A756, 419 (2005) [arXiv:hep-ph/0411405]; Phys. Lett. B610, 253 
(2005) [arXiv:hep-ph/0501193]. 

[21] A. Kovner and M. Lublinsky, Phys. Rev. D 71, 085004 (2005) [arXiv:hep-ph/0501198]. 

[22] Y. Hatta, E. Iancu, L. McLerran, A. Stasto and D. N. Triantafyllopoulos, Nucl. Phys. A 764, 423 (2006) 
[arXiv:hep-ph/0504182]. 

[23] E. A. Kuraev, L. N. Lipatov, and F. S. Fadin, Sov. Phys. JETP 45, 199 (1977); Ya. Ya. Balitsky and L. N. Lipatov, 
Sov. J. Nucl. Phys. 28, 22 (1978). 

[24] J. Bartels, M. Braun and G. P. Vacca, Eur. Phys. J. C40, 419 (2005) [arXiv:hep-ph/0412218] ; J. Bartels and 

C. Ewerz, JEEP 9909, 026 (1999) [arXiv:hep-ph/9908454] ; J. Bartels and M. Wusthoff, Z. Phys. C66, 157 (1995) ; 
A. H. Mueller and B. Patel, Nucl. Phys. B425, 471 (1994) [arXiv:hep-ph/9403256]; J. Bartels, Z. Phys. C60, 471 
(1993). 

[25] M. A. Braun, Phys. Lett. B632 (2006) 297 [arXiv:hep-ph/0512057]; arXiv:hep-ph/0504002 ; Eur. Phys. J. C16, 337 
(2000) [arXiv:hep-ph/0001268]; Phys. Lett. B 483 (2000) 115 [arXiv:hep-ph/0003004]; Eur. Phys. J. C 33 (2004) 113 
[arXiv:hep-ph/0309293]; Eur. Phys. J. C6, 321 (1999) [arXiv:hep-ph/9706373] ; M. A. Braun and G. P. Vacca, Eur. 
Phys. J. C6, 147 (1999) [arXiv:hep-ph/9711486]. 

[26] E. Iancu, A. H. Mueller and S. Munier, Phys. Lett. B 606 (2005) 342 [arXiv:hep-ph/0410018]; E. Brunet, B. Derrida, 
A. H. Mueller and S. Munier, arXiv:cond-mat/0603160; Phys. Rev. E73 (2006) 056126 [arXiv:cond- mat/0512021]. 

[27] R. Enberg, K. Golec-Biernat and S. Munier, Phys. Rev. D72 (2005) 074021 [arXiv:hep-ph/0505101]. S. Munier, 
arXiv:hep-ph/0608036. 

[28] M. Kozlov, E. Levin and A. Prygarin, arXiv:hcp-ph/0606260. 

[29] P. Grassberger and K. Sundermeyer, Phys. Lett. B77 (1978) 220. 

[30] K. G. Boreskov, "Probabilistic model of Reggeon field theory," arXiv:hep-ph/0112325 and reference therein. 
[31] A. H. Mueller, Nucl. Phys. B415, 373 (1994); ibid B437, 107 (1995). 

[32] V. A. Abramovsky, V. N. Gribov and O. V. Kancheli, Yad. Fiz. 18, 595 (1973) [Sov. J. Nucl. Phys. 18, 308 (1974)]. 

[33] A. H. Mueller and B. Patel, Nucl. Phys. B425, 471 (1994); A. H. Mueller and G. P. Salam, Nucl. Phys. B475, 293 
(1996), [arXiv:hep-ph/9605302]; G. P. Salam, Nucl. Phys. B461, 512 (1996); E. Iancu and A. H. Mueller, Nucl. Phys. 
A730 (2004) 460, 494, [arXiv:hep-ph/0308315],[arXiv:hep-ph/0309276]. 

[34] Y. V. Kovchegov and E. Levin, Nucl. Phys. B 577 (2000) 221 [arXiv:hep-ph/9911523]. 

[35] E. Levin and M. Lublinsky, Nucl. Phys. A730, 191 (2004) [arXiv:hep-ph/0308279] . 

[36] E. Levin and M. Lublinsky, Phys. Lett. B607, 131 (2005) [arXiv:hep-ph/0411121]. 

[37] E. Levin and M. Lublinsky, Phys. Lett. B 521 (2001) 233 [arXiv:hep-ph/0108265]; Eur. Phys. J. C 22 (2002) 647 
[arXiv:hep-ph/0108239]. 

[38] S. Bondarenko, E. Gotsman, E. Levin and U. Maor, Nucl. Phys. A 683 (2001) 649 [arXiv:hep-ph/0001260]. 
[39] Y. V. Kovchegov, Phys. Rev. D 72 (2005) 094009 [arXiv:hep-ph/0508276]. 

[40] K. G. Boreskov, A. B. Kaidalov, V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C 44, 523 (2005) 
[arXiv:hep-ph/0506211]. 

[41] M. Ciafaloni and G. Marchesini, Nucl. Phys. B 105, 113 (1976). 

[42] M. Kozlov and E. Levin, Nucl. Phys. A739 (2004) 291 [arXiv:hep-ph/0401118]. 

[43] I. Gradstein and I. Ryzhik, " Tables of Series, Products, and Integrals", Verlag MIR, Moskau,1981. 

[44] A. Kovner and M. Lublinsky, JHEP 0611 (2006) 083 [arXiv:hep-ph/0609227]; A. Kovner, M. Lublinsky and 
H. Weigert, arXiv:hep-ph/0608258 J. Jalilian-Marian and Y. V. Kovchegov, Phys. Rev. D 70, 114017 (2004) 
[Erratum-ibid. D 71, 079901 (2005)] [arXiv:hep-ph/0405266]. 



- 20 - 



